Singularity-free Cosmo-logical Solutions with Non-rotating Perfect Fluids 
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The paper establishes the result that solutions of the type described in the title of the article are 
only those that have been already presented in the literature. The procedure adopted in the paper 
is somewhat novel - while the usual practice is to display an exact solution and then to examine 
whether it is singularity free, the present paper discovers the conditions which a singularity free 
solution of the desired type must satistfy. There is no attempt to obtain exact solutions. Simply, 
the conditions that were ad-hoc introduced in the deduction of singularity free solutions are here 
shown to follow from the requirement of non-singularity. 
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I. INTRODUCTION 

A number of singularity theorems proved in the late 
sixties of the last century and afterwards led to a belief 
that a singularity free cosmological solution in general 
relativity is not possible unless one violates the strong 
energy condition and/or causality. Here we are not in- 
terercstcd in a critique of these theorems or of the justifi- 
cation of the belief about the non-existence of singularity 
free solutions. We simply note that when in 1990 Sen- 
ovilla |Q presented a non-singular solution with a perfect 
fluid having p = 3p > 0, it appeared to be something 
exceptional. Shortly afterwards, Ruiz and Senovilla 0] 
presented a family of non-singular solutions for perfect 
fluids with or without an equation of state. None of 
these investigations were motivated by a specific search 
for non-singular solutions but came as a by-product of 
an exhaustive investiagtion of G2 cosmologies with the 
following mathematically simplifying assumptions: 

(a) The two killing vectors commute and are 
hyper-surface orthogonal. 

(b) The metric tensor components are separa- 
ble functions i.e., one can write g^ v — T^R^ 
where T's and R'S are functions of the time 
coordinate t and the space coordinate r re- 
spectively. The coordinates t and r are re- 
spectively the lines tangential to the velocity 
vector and the space coordinate orthogonal 
to the spaces spanned by the killing vectors. 

(c) All the functions T M „ are expressible as 
powers of a single function of t. 

Ruiz and Senovilla did not assume cylindrical symme- 
try but noted with surprise that all singularity free so- 
lutions discovered by them had this symmetry and won- 
dered whether "this type of symmetry could have some 
relevance to the avoidance of singularity." 

Somewhat later, Mars || gave up the assumption that 
the killing vectors are hyper-surface orthogonal and used 
an " ansatz of separation of variables in comoving coordi- 
nates" analogous to the assumption (b) above. However, 
the condition of separatibilty as usually understood does 



not fully hold good in his solution. [The exception is 
.900 = r 2 cosh(2at) + a 2 r 4 cosh - (2at).] Mars's solution 
had the equation of state p — p. Ruiz and Senovilla have 
clearly demonstrated that the only non-singular solutions 
satisfying their assumptions are given by their equation 
(58), Mars's work showed that the generalization to non- 
orthogonal G2 adds merely one member to the family of 
non-singular solutions with non-rotating perfect fluids. 

It is not unreasonable to think that these non-singular 
solutions (discovered, so to say accidentally) are just a 
subset of measure zero amongst an infinite family of non- 
singular solutions - the discovery of the particular ones 
has been facilitated by their simple mathematical form 
originating from the very strong assumptions of Ruiz and 
Senovilla. Such a view would have been confirmed if dur- 
ing the last decade some non-singular solution outside the 
already discovered family were found. But as nothing like 
that has happened, it is tempting to nurse the alterna- 
tive idea that the already discovered solutions exhaust 
all possible non-singular solutions with non-rotating per- 
fect fluids irrespective of any other consideration. The 
present investigation was motivated by such an idea and 
establishes this result by deducing all the assumptions of 
Ruiz and Senovilla as well as Mars from the condition of 
non-singularity. But first let us spell out the meaning of 
different terms used by us: 

(i) Non-singularity: It is meant to indicate that none of 
the physical or geometrical scalars blow up even at infin- 
ity. The metric tensor components may blow up (indeed 
some do blow up in Ruiz-Senovilla and Mars solutions) 
at spatial or temporal infinity but they must be bounded 
and continuous in the finite region of space-time. We 
make no use of the sophisticated idea of geodesic com- 
pleteness in consideration of singularity. Rather we feel 
that identifying singularities with incompleteness of null 
and/or timelike geodesies is in the ultimate analysis 'a 
piece of opportunism as it allows many theorems to be 
proved' 0. 

(ii) Non-rotating: The fluid flow is assumed to be hy- 
drodynamic so that one can introduce a comoving co- 
ordinate system with the velocity vector hyper-surfacc 
orthogonal, i.e., w M = t^ 010 v v v a -^ = 0. 
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(iii) Perfect fluid - The energy stress tensor is of the 
form T£ = (p + p)v^v v — p8^ in the usual notation. This 
means that we are just leaving out viscosity, heat flow 
etc so that the solution presented by Dadhich et al ||] 
are outside our purview. 

(iv) Cosmological: The energy tensor obeys the strong 
energy condition - in particular for the perfect fluid we 
take somewhat stronger constraint < p < p. The pres- 
sure and the density must not vanish anywhere except at 
infinity. We thus exclude the case of a complete vacuum 
or a bounded distribution of matter. We shall exclude 
any discontinuity in T M „ - although that is not inconsis- 
tent with absence of singularities, such discontinuities are 
unusual for simple cosmological models. This will allow 
us to use a single coordinate system as globally valid. 

(v) We further introduce the physical assumption that 
changes of p and p are in the same direction for any 
change of the coordinates. Such a condition seems natu- 
ral for any physically acceptable model. 

It is easy to see that in the absence of rotation, a non- 
singular cosmological solution of the desired type must 
have both acceleration and expansion non- vanishing. 
For, in the absence of acceleration, expansion will nec- 
essarily be associated with a collapse singularity as in 
the Friedmann models and the expansion free case is not 
of interest in cosmology. (They are usually studied as 
stellar models.) 

Dadhich |6) has claimed that instead of all the assump- 
tions of Ruiz and Senovilla, if one, besides our assump- 
tions, assumes that in the comoving coordinate system, 
the metric tensor is diagonal and separable into functions 
of time and space coordinates, one is uniquely led to the 
family discovered by Ruiz and Senovilla. But the addi- 
tional assumptions of Dadhich are quite strong and leave 
out Mars's non-diagonal metric. Our treatment is free of 
these shortcomings and assumes nothing besides the five 
conditions stated above. 

It may be of some interest to note a paper by Senovilla 
0. The Senovilla paper adopts the metric: 

dg 2 = T 2(l+„) E 2n(n-l) ( _ dt 2 + ^2) + ^(l+n)^^^ 
+T 2(l-n)^2(l-n) dz 2 

where T is a function of t alone and S is a function of 
r alone. This metric satisfies all the conditions used by 
Ruiz-Senovilla in their deduction of the singularity free 
perfect fluid solutions, which under the circumstances are 
shown to be only such solutions. Specifically the assump- 
tions are existence of G2, orthogonal separable metric, all 
time functions in the metric are powers of a single func- 
tion. Thus any non-singular solution (they are not ex- 
plicitly displayed) that |7| may imply which differs from 
that of the Ruiz-Senovilla will either violate the perfect 
fluid condition (note that Senovilla allows p r 7^ p z ) or 
violates the strong energy condition. (Senovilla himself 
points out that there are "almost FLRW model" which 
are singularity free but "violate energy conditions".) 



II. THE COORDINATE SYSTEM AND THE 
FUNCTIONS P, p, V„ AND 6 

The non-rotating condition allows us to write the met- 
ric in the form: 



ds 2 



g 00 dt 2 + g ik dx l dx k (i, k =1,2, 3) 



with the fluid velocity components 

1 



v* = 0, v° = 



The acceleration vector has the components 
«« = -^(ln5oo),« ;v = 0. 



(1) 



(2) 



(3) 



Thus the acceleration vector is confined to the three space 
orthogonal to w M and in that space it is a gradient vector. 
Hence we can align a space-like coordinate (say x 1 = r) 
along this vector which will be orthogonal to the other 
three coordinate lines and hence (TO) will be reduced to 



ds 2 = goadt 2 + g\\dr 2 + g a bdx a dx b 



with 



0, 



2,3) 



(4) 



(5) 



We shall consider the domain of the coordinates to ex- 
tend from —00 to +00. In particular we note that as the 
tangent to the r coordinate lines is a gradient vector, any 
r-line cannot form a closed loop. Thus r has to run from 
— 00 to +00 (or in case the r-lines diverge from a point 
as the radial lines in spherical or cylindrical symmetry, r 
will run from zero to 00. 

One may wonder whether the metric form (^) is glob- 
ally valid with a single coordinate system. However, we 
note that our condition (iv) setting our idea of a cosmo- 
logical solution ensures this global validity. Obviously, 
because of (0) and (||), 1700 is a function of r and t only 
(not of x 2 ,jr). From the relation — 0, we get for 
perfect fluids, 



P,i 



p + p 



(In 5oo), i 



(6) 



Hence p is also a function of r and t only. Writing equa- 
tion (H) in the form 



(P + P) = - 



2Pn 



(In 500), 



one sees that p is also a function of r and t only. From 
the other divergence relation 



1 



-P 



P + P 



(7) 



one sees that 9 is also a function of r and t only. Summing 
up, the scalars p, p, 9 and the metric tensor component 
goo (whose gradient determines the acceleration) are all 
functions of r and t only. 
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III. THE VANISHING OF SOME FIRST 
DERIVATIVES 

In two communications , the present author estab- 
lished that for non-singular, non-rotating cosmological 
solution, both the space-time average and the space av- 
erage (defined in a suitable manner), of each of the scalars 
that appear in the Raychaudhuri equation vanish. These 
results in turn require that for distributions of perfect 
fluid which are not bounded, these scalars must vanish 
both at spatial and temporal infinity. Thus p, p, a, 9, v^;^ 
all vanish at r, t — » ±00. Consequently the positive defi- 
nite quantities likep, p, a 2 , 9 2 must have at least one max- 
imum both in r and t. 

Consider now the case of pressure p. As it is a function 



of r and t, the maximum in p in r defined by 



Or 



and 



0$ < will be a line in the r, t space. Along this line we 
shall have, 



d 2 p 
dr 2 



drj 



d 2 p 
drdt 



dt 1 



(8) 



The subscript L indicating that differentials along the 
line are to be taken. Also as ^3 < all along the line, 
the integral J i^dr^ will diverge if the domain of in- 
tegration be infinite but J -g^drr, on integration gives 
(ff )l and will be finite for all limits of integration. Thus 

the only way to satisfy (H|) is to have drj, = and 
also zero at all points on the line. Thus this line is a 
i-line defined by a constant value of r. 

An exactly similar analysis leads to the conclusion that 
(^|) = line is a i-constant line, i.e., a r coordinate line 

and over it also jff§i vanishes, choosing the intersection 
of these two orthogonal lines as the origin (r = 0, t = 0) 
we may sum up as follows: 



dp 
dr 



for r = 



for t = 



(9) 



of p. However, Such a situation will go against our as- 
sumption (v). Hence we have ro = to = 0. 

Next consider the case of (In goo)- In view of equation 
(|), the zero of §f is also a zero of d % 9 r oo) . 

From (H), we get, 



dp/dr 



P 



(1 



p g(lngoo) 
p dr 



If the integral extends from — 00 to +00, p 
. p 



as (1 + is finite and positive everywhere, 500 



0, and hence 
-00 as 

r — > ±00. Again, from (0), by a suitable similar integra- 
tion, \det gik\ = \h\ — > +00 as t — > ±00. But 9 and v^;^ 
vanish at the infinities, consequently goo ~~ * +°° also at 
t — * ±00 and J h | — ► 00 as r — > ±00. Hence both \h\ and 
goo must have minimum with respect to both r and t. 

That these minina will occur at r = and t — again 
follows from (g), (0), (g) and jlO|). A further investiga- 
tion of the vanishing of cr 2 and u^;^ at the infinities show 
that g rr and 7 (the two dimensional metric determinant 

= t^t) also blow up at the infinities. 

laid' ^ 

Thus we get 



2{;J u _ d(lngoo) 



dr 



= for r = 



d (In goo) 



(9/- 



= for t = 



(11) 



d 2 (In goo) 



drdt 

Similarly from (^), we get, 
d1n\h 



= if either r = or t = 



9r 



= for r = 



(9t 



(12) 



9 2 p 



drdt 
Exactly similarly. 



= if either r = or t = 



dp 

— = for r = constant = ro (say) 
dr 



dp 

dt 



for t = constant = to (say) 



(10) 



d 2 p 
drdt 



= if either r = ro or t = to 



d 2 In ft 
drdt 



= if either r = 



t = 



The pair of orthogonal lines r — and t — are the 
lines of vanishing of first derivatives of p, p, In goo, In ft. It 
is easy to see that at no other value of r, can ^ vanish. If 
possible let us suppose that at r — 0), |^ = 0. From 
our previous analysis, it is clear that non-singularity re- 
quires that this cannot be an isolated point but would 
be the entire t coordinate line at which r = R. Also 
as r = R follows the maximum of p at r — 0, will 
benegative or zero over the entire line. Consider now the 
Raychaudhuri equation, 



We have introduced ro , to to indicte the possibility that 
the maximum of p may not coincide with the maximum 



2a 2 + ±9 2 + ^-(p + 3p) = v»- 1 ,~t 



(13) 
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where 



and 



1 88 

/ffoo dt 



Now all the terms of the left hand side of (|13| ) are pos- 
itive and the first term on the right is zero or negative as 
we have just seen. Hence 9 is negative at all points on 
the line and 
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In 



dt>0 



which contradicts the condition of regularity that 9 is 
zero at both t — > ±oo. 

In a similar manner it is easy to see that there is no 
other zero of ^ except at t = 0. 



IV. EXPRESSIONS OF THE DIFFERENT 
VARIABLES IN TERMS OF SEPARATED 
FUNCTIONS 

The results of the last section allow us to express p in 
the following form 

p = ER i T i Q(r, t) 

where R's are linearly independent functions of r having 
= at r = and vanishing at r — » ±oo. T's are simi- 
larly functions of t alone. We have put in the summation 
to take care of the possibility that there may be more 
than one function satisfying the imposed conditions. (In- 
deed if f(r) is one such function then any positive power 
of f(r) will satisfy the same condition.) Q(r,t) is any 
function involving both r and t (if possible). However, 
as there is no other zero of ^ and except those cov- 
ered in R's and T's, the drivarive of Q(r, t) with respect 
to r and t cannot have any zero. Cosequently, if Q is 
to be bounded everywhere, it must be a constant and 
can be absorbed in the R's and T's. (It is assumed that 
the summation may be over infinite terms to include all 
possible R's and T's.) 

The functions R and T will be determined by the 
field equations which are second order partial differen- 
tial equations - hence there can be at most two linearly 
independent R functions and two T functions. So finally, 



p = RpTp 



and similarly 



RpTp 



RpTp 



RpTp 



m 3oo = RqTq + RqT 



\nh = R h T h + R h f h + M(x a ) 



(14) 



(146) 



(14c) 



(14d) 



where the unbarred and barred functions with the same 
subscript are linearly independent. 

As yet, we have not proved that h is independent of the 
coordinates x 2 , x 3 , hence we have added the term M(x a ) 
in (14d). No such term is required in the other three ex- 
pressions as p, p,ln<7oo have been shown to be functions 
of r and t only. 

Putting in the expressions (14a, b,c) in equation (6), 
we get, 



R'pTp+R'pTp — —-(RpTp+R p Tp+RpTp+RpTp)(R' To+R' T ) 

(15) 

where superscript primes indicate differentiation with 
respect to r. 

Using now the non-linear dependence of the barred and 
unbarred functions, we get from equation (15), 



In 500 = Rq + T 
p = aRpTp + bRpTp 
P — RpTp ~\~ RpTp 



(16a) 
(166) 
(16c) 



where a and 6 are constants and Rq , T are different 
from those in (14c). Also, 

In R v = — — ^— In R n 



and 



1 + a 



Rn — In Rn — -lnR n 



\ ■ a 1 • 6 ■ 

Similarly plugging in (14d) in equation (7), we get, 



ln\h\ = R h + T h + M(x a ) 
1 + ba 



lnT n 



1 — a 6 



In XL 



1 + a 



(17a) 
(176) 

(17c) 



We have already noted that gn and 7 also blow up at 
the infinities of r and t. Hence they also must be express- 
ible in form similar to the above. Keeping in mind the 
relation \h\ = |7||<?n|, we can write, 



ln|ffii| =R 1 +T 1 + N{x a ) 



(18a) 
(186). 



ln| 7 | =i? 7 + T 7 + P(a; Q ) 

Note that we can make transformations t — fit) and 
f = 4>(r) to have 

ln|sm| - In .goo = N(x a ) 

or, 3n/(?oo = a function of x a . 
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V. PROOF OF THE EXISTENCE OF G 2 



9ll,a = 



We shall proceed in several steps. The first step is to 
prove that in the solutions that we are seeking, shear 
must be non-vanishing. 

For this it is necessary to use the field equation 



R 0l = -8tt(T - T 0l ) 



(19) 



The above equation for a perfect fluid in irrotational 
motion can be written in the form 



In our present coordinate system this becomes 



i\k 



(20) 



(21) 



where \k indicates covariant differentiation with the three 
space metric. If now a\ vanishes, 9 will be spatially uni- 
form and then either the situation is static (6 — 0) or the 
solution will have singularities. 

If the three non-zero values of shear eigenvalues are 
Si, 5*2 and S3, we get 

^m^-km^ h ^ (22) 

i.e., it is a function of r and t only. 

In the next step we show that a 2 and u;^ are inde- 
pendent of the coordinates x a . Looking back at equation 
(13), we see that in view of (16) and (17) all the terms 
except 2tr 2 and v;j[J are functions of r and t only. Hence 
we must have, 



But 



(2a 2 ), a =(v;^), a 



■-(vi^V-ij),! 



,N/2 



9Q0 



-[(In 3oo) 



1 



(In 5oo), 1 (hi7),i 



Excepting e N / 2 , g 00 and the term within the braces are 
functions of r and t only and these have no zero for finite 
values of r and t. Thus, 



We now come to the last step of our proof of existence 
of G 2 - For this we introduce three orthonormal vectors 
\\\, XL, A31, the three unit eigenvectors of the shear ten- 
sor. Using the following formulas given in Eisenhart's 
book @ 

^;fc A l| A f| A m| = (S h - Si) lhVm (h I) 



r ■ ,y \ 3 \ k - \ kdSh 



we get from (21), 



[Si,i + tt-(5i - 5 2 )7i22 - 



[5*2.2 + Tt(S2 
A |2 



1 



#3)7133] (24) 



5*3)7233 + to- (S2 - S'i)7 2 ii] (25) 



A |2 



= [5*3,3 + TT( S 3 - 52)7322 + Tt(53 - 5i)7; 



A |3 



7311J 



(26) 



where 7's are Ricci relation coefficients. As <7n, a = 0, A 
is geodetic in the three space. Hence 7211 = 7311 = I 
Now from (25) and (26) as S ha = we get, if 5 2 ^ 5 3 



7233 = 7322 = 



(27) 



If 52 =5*3, then the existence of G*2 can be readily proved 
but as will be apparent from our later discussion, this 
case will lead to the vanishing of shear and hence to sin- 
gular solutions. From (27), owing to the antisymmetry 
labc = —Ibac, it follows that all the 7's involving the sub- 
scripts 2 and 3 vanish and consequently the two spaces 
spanned by x 2 and x 3 are all flat; i.e., we may choose 
coordinates such that g a b,c = 0. Thus all the metric ten- 
sor components are independent of x a - the space time 
admits translations along x 2 and a; 3 . 



VI. CYLINDRICAL SYMMETRY 



From the field equations, one can deduce in case of 



(2a 2 



N/2 . . . irrotational motion of a perfect fluid (cf. Ellis [fLl| 

= (e ' ), a x a non — vanishing expression depending on r and t. 4 — ' 



will 



Consequently, if TV is not a constant (2cr 2 ), Q or (u;^ . 
have no zero in the finite region of r and t. But as a 2 
vanishes at t = for all values of r and x a and (2<r 2 ), a 
or ({i;^ , a also vanish there, and it follows that N is a 
constant and hence 

5i, a = 5 2>a = 53, a = (23) 



-9 2 + 2<t 2 + mwp 



(28) 



where *R is the Ricci scalar for the three space t 
Now by direct calculation, 



const. 



1 8e 1 n n 

V~9n or 2 



(29) 
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where 9 and E are the 'expansion' and 'shear scalar' for 
the unit space-like vector £ l : 



—Si, i.e., 



© = a 



i 



V-9r- 



ny2 v~i y^t/c 

ZZj — t-'ik 2 -' 



^ik - ^(^k + (6>|i) - ^(^ifc + teffc)^ 

where |i indicates covariant differentiation with the three 
spaces metric. The reader may note the similarity be- 
tween equation (29) and the Raychaudhuri equation in 
the 4-dimensional space time. Again with G2, 



"R =* R\ 



1 



de 



y/-gil dr 



e 2 



(30) 



Combining (28), (29) and (30), one gets, 



— -?e 2 = -V + 2a 2 + 16np + 2E 2 (31) 
V-gii or 2 3 

For the three space at t = 0, 9 2 = a 2 = and we get 
from (31) 



99 



-9 2 ) > 



(32) 



The inequality (32) shows that 9 would blow up leading 
to the collapse of the two-space (i.e., 7 = 0). This will 
not mean a singularity iff x 2 or x 3 or both be angle coor- 
dinate so that 7 will naturally vanish at r = 0. However 
both cannot be angle coordinates for then E 2 would van- 
ish at r = and (31) would make p blow up at that point. 
With one only of x 2 , a: 3 being angle coordinate both the 
left hand side as well as E 2 would blow up leaving a finite 
difference for p. 

Thus one of the translations is along an angle coordi- 
nate indicating that we have cylindrical symmetry. 



VII. RELATION BETWEEN DIFFERENT TIME 
FUNCTIONS 

Let us identify x 2 as the angle coordinate. Then the 
condition of elementary flatness gives S\ = S2. Plugging 
in this result in (24), we get 



or 



Si,i + 3y(lng 33 ),i = 



1 Si 1 

-5i(lnpoo),i -3y (hic/ 33 ),i = -6»(ln.9 00 ),i 



so that 



gngggk _2 j^-f) 



(In. g o), 1 3 Si 



(33) 



From what we have proved it is clear that the left hand 
side of (33) is a function of r alone while the right hand 
side a function of t only. Hence both sides must be con- 
stant. So that we have 



£2 = Si = nd; S 3 



-2nd 



(34) 



where n is a constant. Equation (34) leads to the de- 
sired result that Ti = T2 = Tq, and Th are expressible 
as powers of a single time function. Equations (17b) and 
(17c) then show that T p and T p are also of the same type. 
Hence we have deduced the assumption (c) of Ruiz and 
Senovilla. 

At this stage it is clear that we have been able to de- 
duce all the assumptions made by Ruiz and Senovilla 
except the hypersurface orthonormality of the killing vec- 
tors. This is as it should be for Mars has shown that non- 
singular solutions do exist where the killing vectors are 
not hypersurface orthogonal. Also, our demonstration of 
separatibility stops at that of 7 and does not include the 
non-diagonal g a i,. This again is consistent with Mars's 
finding. 

Thus we may claim to have achieved our objective 
namely that with our conditions (i) to (v), the already 
known solutions given by Ruiz-Senovilla and Mars are 
the only possible non-singular solutions. 



VIII. CONCLUDING REMARKS 



It is tempting to ask whether the present work can 
be extended to more general cases, i.e., whether one can 
establish specific conditions about non-singular cosmo- 
logical solutions when one or more of the conditions of 
irrotational motion, perfect fluidity and |^ > is given 
up. Already there exist quite a number of imperfect fluid 
non-singular solutions but they vary widely in their char- 
acteristics. Thus besides cylindrically symmetric solu- 
tions, there are also spherically symmetric solutions |l^] 
and even an oscillatory solution pi . In fact for imperfect 
fluids, the divergence relations = written out ex- 
plicitly are much more complicated than (6) and (7) and 
do not seem amenable for simple general analysis, again 
if rotation be present, non-diagonal terms gm appear in 
the metric tensor and further calculation becomes much 
involved. So far as the condition ^ > is concerned, the 
present author believes that it can be done away with al- 
though right now he is not in a position to show that. 
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